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Abstract 

Cosmic rays (CR), constrained by scattering on magnetic irregularities, are be¬ 
lieved to propagate diffusively. But a well-known defect of diffusive approximation 
whereby some of the particles propagate unrealistically fast has directed interest to¬ 
wards an alternative CR transport model based on the “telegraph” equation. However, 
its derivations often lack rigor and transparency leading to inconsistent results. 

We apply the classic Chapman-Enskog method to the CR transport problem. We 
show that no “telegraph” (second order time derivative) term emerges in any order of 
a proper asymptotic expansion with systematically eliminated short time scales. Nev¬ 
ertheless, this term may formally be converted from the fourth order hyper-diffusive 
term of the expansion. But, both the telegraph and hyperdiffusive terms may only 
be important for a short relaxation period associated with either strong pitch-angle 
anisotropy or spatial inhomogeneity of the initial CR distribution. Beyond this period 
the system evolves diffusively in both cases. The term conversion, that makes the tele¬ 
graph and Chapman-Enskog approaches reasonably equivalent, is possible only after 
this relaxation period. During this period, the telegraph solution is argued to be un¬ 
physical. Unlike the hyperdiffusion correction, it is not uniformly valid and introduces 
implausible singular components to the solution. These dominate the solution during 
the relaxation period. As they are shown not to be inherent in the underlying scattering 
problem, we argue that the telegraph term is involuntarily acquired in an asymptotic 
reduction of the problem. 


1. Preliminary Considerations 


The problem addressed here is fundamental and not new to the cosmic ray (CR) transport 
studies. It can be formulated very plainly: How to describe CR transport by only their isotropic 
component, after the anisotropic one has decayed by scattering on magnetic irregularities? 
Suppose the angular distribution of CRs is given by the function f(/j.,t,z) obeying an equation 
from which the rapi d gyro-p hase rotation is already removed (drift approximation, e.g., 
Yedenov et al. 19621: Kulsrud 2005 ) 




d£ 

dfl' 


( 1 ) 


Here z is the local coordinate along the ambient magnetic field, /i is the cosine of the particle pitch 
angle, and is the pitch angle diffusion coefficient. Now, we make the next step in simplifying 
the transport description and seek an equation for the pitch-angle averaged distribution 
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fo(t,z) = - / f(n,t,z)dn = {f). 


The basic solution to this problem has been known for at least half a century (e.g., Jokipii 1966 
and references therein). To the leading order ml/3$ (assuming the characteristic scale and time 
of the problem being longer than particle mean free path and collision time) it can be obtained 
straightforwardly by averaging cq.([Q) 


dfo 

dt 




and substituting df/dfl <C fo from eq.© as: 


d£ 

dfl 


v df 0 


23$ dz 

Thus, the following diffusion equation results for fo : 


( 2 ) 


dfo 


d / 1 -M 2 \ dfo 


dt 4 dz 


/ dz ' 


(3) 


A questionable point of course is neglecting df/dt in favor for vdf/dz in eq.©. It is 
somewhat justified by the small parameter & 1 <C I in the, final result, given by eq.©, making 
df/dt hopefully small. On the other hand this is true for dfo/dt but not necessarily for df/dt, 
since the latter may contain also the rapidly decaying anisotropic part f — f — fo of the initial CR 
distribution. For 3ft > 1, however, / must die out and neglecting df/dt appears plausible for the 
long-term CR transport. At the same time, dfo/dt is large when fo is very narrow in z initially, 
such as in the fundamental solution. In the sequel, these aspects of the CR propagation will be a 
key for choosing an appropriate asymptotic reduction method. 

However convincing the justification, the CR diffusion model encounters the problem of a 
superluminal, or simply “too-fast” particle propagation. Although rather common for diffusive 
models, the problem is largely ignorable as long as the number of such particles remains small. 
There are cases, however, such as the pro pagat ion of ultra high-energy cosmic rays, where this 
problem mus t be addressed ( Aloisio et al. 2009h . Various attempts, starting as early as in 60s, e.g., 
( Axfordl 1965 ). have been made to devise a better transport equation for CRs. Unfortunately, in 
our view, they lack mathematical rigor and clarity and sometimes lead to inconsistent results. 

In the most recent telegraph model, due to Litvinenko & Schlickeiserl ( 2013 1. a higher order 
in \/3$ -C 1 term was included by retaining d f /dt, dropped in the simplest derivation above. This 
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strategy gave rise to an additional d 2 fo/dt 2 -term in the “master” equation. This additional term 
transforms eq.© into a “telegraph” equation: 


with 


dfo , ^ d 2 fo d df 0 k df 0 


dt 


dt 2 dz dz L dz 


(4) 


v 2 /l-At 2 

k= j 


T = 


d\x/3> 



, ZT 1 = -B~ l dB/dz 


(5) 


For the sake of comparison with earlier telegraph equation results that will be made in Sec HI 
we have added here the magnetic focusing effect (the last term on the r.h.s. with B(z) being 
the magnetic field), not included initially in cq. ([]]). Eq.© is just a linear equation that 
can be solved immediately. The fundamental solution to cq.© that starts off from a S (z) 


distribution, instantaneously released at t — 0, is as f ollows (e.g.. lGoldsteinlll95ll: lAxfordlll965 
Schwadron & Gombosill 1994c Litvinenko et al.ll2015l. L 1 =0, for simplicity) 



Here /o.i denote the modified Bessel functions and H - the Heaviside unit function. 

One promising aspect of the telegraph equation is that it allows for a ballistic mode of CR 
propagation when the initial conditions empower the higher-order derivative terms to dominate 
(at least in early phase of evolution). If, in addition, T has a proper value, the bulk speed of 
CRs may also be realistic. For example, this speed was derived in (Earl 1973 ) to be v/\/3, 


which has also been used earlier by Axford 1965. This is just the rms velocity projection 


of an isotropic, one-sided CR distribution on z-axis, which appears to resolve the issue with 
the superluminal propagation. What is worrisome here is that this bulk speed essentially 
requires a one-sided “isotropic” CR distribution which, of course, is highly anisotropic overall, 
contrary to the basic assumption of most treatments. So, we need to start with isotropic initial 
distribution but, taking it narrow and symmetric in z (say Gaussian), the higher-order derivative 
terms will, again, dominate in eq.@ and the single CR pulse will split into two, propagating 
in opposite directions at the speeds ±y/k/T (as dfo/dt — 0 at t — 0, according to eq.lfll). 
Neglecting dfo/dt (which is justified for /o sufficiently narrow in z), the solution is simply 
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fo(z,t) = F (z— y/k/Tt j + F (z + a Jk/Tt'j , where IF (z) — f(z , 0). This result casts doubts 
on whether the telegraph term can be dominant under the assumption of frequent CR scattering 
(asymptotic expansion in small 1 /&), as the initially sharp profile does not spread. Bringing 
df 0 / dt back into the equation will only damp but not spread the profile, as clearly seen from the 

solution given in eq.©, where F (z ± y/k/Tt ) = 8 (z ± y/k/Tt ). Besides, the bulk CR speed 


derivation of the telegraph e quation for an arbitrar 


is not universally accepted. iGombosi et ajJ (1 19931) : iPauls et al. 


V 3>{u) by 

Litvinenko & Schlickeiser 

(2013 

), 

Pauls et al. 

(1993 

) and 

Schwadron & Gombosi 


( 19941) . using simplified forms of D (p), advocate the value y/5j\lv for the propagation speed. 


The last result is consistent with our calculations below, but with strong reservations 
regarding the telegraph equation set out later in the paper. Here we merely note that the solution of 
telegraph equation specifically considered by Litvinenko & Schlickciscrl ( 2013 ). which does not 
have the property of splitting the initially narrow pulse into two, does not conserve the total CR 
number N — J fodz- It starts off from N — 0 which is unphysical, as the equation has no source 
on its r.h.s. To conserve N, two <5— pulses in eq.© are necessary. Those have been added to the 
treatment by Litvi ne nko & Noble 120131) . but the 5— pulses have not been shown on their plots, 
for obvious reason. So, the comparison of this solution with the solution of the original scattering 
problem is rather misleading. The disagreement on the propagation speed y/k/T is also critical 
as the solution in eq.© is cut off at a point z moving with this speed. For z < y/k/Tt the profile 
is close to a Gaussian (for t T, where T is the scattering time), so small variations in the speed 
can produce significant variations in the solution. We will also return to this later. 

Another disadvantage of the telegraph equation © is that it is no longer an evolution equation 
and requires the time derivative d t fo as an initial condition. Although this can be inferred from 
the angular distribution at t — 0 using eq.©, the “telegraph” description of CR transport is not 
self-contained. We show below that the T -term in eq.® is subdominant in an asymptotic series 
for Q)t > 1, thus representing transients in the CR transport. Strictly speaking, it should be omitted 
in the asymptotic transport description along with the small hyper-diffusion term S*h/3z\ 
particularly if the term ~ d 3 fo/dz 3 does not vanish. The latter was not included in eq.©, as it was 
obtained by applying insufficient direct iteration to eq.© in dLitvinenko & Schlickeiser 120131) . 
or assuming symmetric scatterin g, S’ (-u) — St (jl ), when it indeed vanishes. This symmetry 
restriction was relaxed in (Pauls et al. 1993). 


The reasons why we undertake the derivation of master equation to a higher (fourth) order 
of approximation for an arbitrary @>(n) are several. First, it is necessary to clarify the role of 
the telegraph term entertained in the literature as an allegedly viable alternative to the standard 
diffusion model. Second, it is important to obtain the transport coefficients valid for arbitrary 
f^(ju), that is for an arbitrary spectrum of magnetic fluctuations. As we will show, the previous 
such derivation due to lLitvinenko & Schlickeiser] 120131. does not include the third order term, 
while including only one fourth order term, while there are more such terms (see eq. If27ll i. 
Furthermore, the diffusion equation © supplemented by a convective term u (z) dfo/dz for the 
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case of the bulk fluid (scattering center) motion with velocity u, has long been and remains the 
main tool for diffusive shock acceleration (DSA) models. An accurate assessment of the next 
non-vanishing term, not included in eq.©, is thus utterly important for the DSA, particularly as 
claims are being made about the necessity to include the telegraph term in the CR transport. In 
most DSA applications, it is crucial to allow not only for an arbitrary fluctuation spectrum *2) but 
for its dependence upon /o as well. This dependence directly affects the particle spectrum and 
acceleration time. We will discuss these aspects briefly in Secj6l 


In the next section, the basic transport equation with magnetic focusing is introduced and the 
shortcomings of a reduction scheme based on direct iterations are demonstrated. The appropriate 
asymptotic method is elaborated in Sec[3] Apart from what we already discussed regarding the 
telegraph equation. , the objective of Sec 0] is to create a frame work suitable also for nonlinear (e.g ., 
Ptuskin et al. 2008 : Maikov et al. 2010b) and quasi-linear ( Fujita et al. 201 ll: Maikov et al.l 2013 ) 


versions of CR transport which are important for both the DSA and for the subsequent escape 
of the accelerated CR. In these settings, the CR pressure is high enough to strongl y modify at 
least the pitch-angle diffusion coefficient 2 and possibly the shock structure itself ( Maiko v et al. 


2010bj) . In Secs I?] and [5] the implications of our results for the telegraph model and for the 


long-time CR propagation are discussed, while Secj6]concludes the paper. 


2. CR Transport Equation and its Asymptotic Reduction 


Energetic particles (e.g., CRs) in a magnetic field, slowly varying on the particle gyro-scale, 


are transported according to the f ollowing gyro-phase averaged equation, e.g. (lYedenov et al. 
1962 : Jokipiil 1966 : Kulsrud 2005 ) 


df df 0 / 2 \ df d / df 

+ +v— (1 -/I 2 ) — = — — vD (/i) (1 -/!-) — 

dt dz 2 v ^ ’ d/j, dpi ^ ’ dll 


(V) 


Here v and /i are the particle velocity and pitch angle, z points in the local field direction, 
a — —B l dB/dz is the magnetic field inverse scale and v is the pitch angle scattering rate, while 
D(li ) 1 depends on the spectrum of magnetic fluctuations. As the fastest transport is assumed 

to be in /i, we introduce the following small parameter 


e 


v A 


( 8 ) 


where A is the particle mean free path and l is a characteristic scale that should be chosen 
depending on the problem considered. One option is the scale of B (z), in which case / ~ 0 ! . 
If the CR source is present on the r.h.s. of eq.flT]]), its scale can be taken as l. Finally, / can be 
the scale of an initial CR distribution. Strictly speaking, the shortest of these scales should be 
taken as /. The problem with the initial CR distribution is that in the most interesting case of the 
fundamental solution this scale is zero. Therefore, over the initial period of CR spreading, before 
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the actual CR scale l (t) ~ // (df/dz) exceeds the m.f.p. A, direct asymptotic expansions in small 
e remain inaccurate. The goal here is to choose the least inaccurate out of all possible expansion 
schemes. At a minimum, it should be the one that does not introduce additional singularities, 
apart from the initial delta function 8 (z), that must spread out under the particle recession and 
collisions. Therefore, while taking / = const in eq.®, and assuming / A, caution will be 
exercised during the initial phase of the CR relaxation when the terms with higher spatial and time 
derivatives are large, even if they contain small factors e' 1 1. By measuring time in v 1 , z in /, 

and simply replacing o7 —>• o ~ 1, the above equation transforms as follows 


d JL 

dt 


d_ 

dfl 




dl 

dfl 



(9) 


A suitable scheme for asymptotic reduction of the above equation using £ <g; 1 is due to 
Chapman and Enskog, suggested in development of the earlier ideas by Hilbert (a good discussion 
of the history of this method with mathematical details is given by Cercignanil 1988). Originally, 
it was applied to Boltzmann equation in a strongly collisional regime. Similar approaches have 
been used in plasma physics^e.g., in regards to the hydrodynamic description of collisional 
magnetized plasmas (IBraginskii 1965b and the problem of run-away electrons ( Gurevichl 1961 : 


Kruskal & Bernsteinl 1964). 


Regardless of the asymptotic scheme, eq.© suggests to seek / as a series in £ 


f — fo + £ /i + e fi + ■ ■ ■ = fo + f 


( 10 ) 


where 


(/Wo, with (•) = ^ f (-)dll, 


( 11 ) 


so that (/) = (/, 7 >o) = 0. The equation for /o, which is the main (“master”) equation of the 
method, takes the following form 


dfo 

dt 


= —£ 


d 

dz + ° 




n= 1 


n ~ l 1 (!-At 2 ) —\ 


dll/ 


( 12 ) 


We see from this equation that, s imilarly to the case of Lorentz’s gas in an electric field (IGurevich 
1961 : Kruskal & Bernstein 1964b . /o depends on the “slow time” = £ 2 / rather than on t. Indeed, 
the two problems are similar in that they describe diffusive expansion of particles in phase space. 
The expansion occurs in z-direction for the CR diffusion problem and in energy for runaway 
electrons. The expansion is driven by a rapid isotropization in pitch angle plus the convection in z- 
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direction, or acceleration in the electric field direction, for the CR transport and electron runaway, 
respectively. 


The slow dependence of /o on time in eq.(fT2l) may suggest to attribute the time derivative 
term in eq.@ to a higher order approximation (thus moving it to the r.h.s.). Such ordering has 
been employed by Litvin enko & S chlickeiser ( 20131) and the term <=c d 2 fo/dt 2 has been produced 
in eq.(fT2l). Obviously, a continuation of this process would result in progressively higher time 
derivatives of /o, corresponding to shorter and shorter times in the initial relaxation. These 
transient phenomena will be removed using the Chapman-Enskog asymptotic reduction scheme 
in the next section. 


Unlike /o, / in eq.tlTOl) does depend on t as on a “fast” time. Therefore, it is illegitimate to 
attribute the first term on the l.h.s of eq.® to any order of approximation different from that of the 
second term, notwithstanding its fast decay for t > 1. Thus, using eqs.OTTOl) we must apply the 
following ordering 


dfn 


dt 


d_ 

<9/i 


DGu)(i-/t 2 ) 




(7 

2 


(i-U) 



(13) 


The above expansion scheme is sufficient to recover the leading order of /o evolution from 
eq.(fl2l) by substituting there <9/i/<9/i — (2D) -1 dfo/dz, obtained from the last equation 

for t > 1. However, this scheme is not suitable for determining /„ for n > 2 to submit to 
eq.(fl2l). Indeed, as it may be seen from eq.(fl3l). the solubility condition for J 2 at / 1 is 

(d/dz + o) ((l -/i 2 ) <9/i/<9/i) « - (d/dz + (j) ((1 -R 2 ) /2D) dfo/dz = 0. This is clearly 
too strong a restriction. The reason for this inconsistency of the direct asymptotic expansion is 
that /o depends on time much slower than f n> 0 , so a slow time (2 = £ 2 t needs to be taken into 
consideration. The Chapman-Enskog method has been developed for such cases, and we will 
make use of it in the next section. 


3. Chapman-Enskog Expansion 


As we have seen, the asymptotic reduction of the original CR propagation problem, given 
by eq.©, to its isotropic part cannot proceed to higher orders of approximation using a simple 
asymptotic series in cq. (171)1) and requires a multi-time asymptotic expansion. In Chapman-Enskog 
method the operator d/dt is expanded instead. Its purpose is to avoid unwanted higher time 
derivatives to appear in higher orders of approximation. This is very similar to, e.g., a secular 
growth in perturbed oscillations of dynamical systems. To eliminate the secular terms, one 
seeks to alter (also expand in small parameter) the frequency of the zero order motion, which 
is similar to the d/dt expansion. One example of such approach may be found in a derivation 
of hydrodynamic eq uations for strongly collisional but magnetized plas mas, starting from 
Boltzm ann equation ( Mikhailovskvl 1967 ). The classical monograph by Chapman & Cowling 
( 199 lh (Ch.VIII) gives another example of a subdivision of d/dt operator for solving the transport 
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problem in a non-uniform gas-mixture. Expanding d/dt operators eliminates secular terms, such 
as the telegraph term. Perhaps more customary today and equivalently is to introduce a hierarchy 
of formally independent time variables (e.g., Navfeh 198 lh t —> to, h,..., so that 


d d d 2 d 

dt dto dt\ dt 2 

Instead of eq.(fl3l). from eq.@ we have 


(14) 


dfn 


dto 


d_ 

dfl 







dfn -1 y df n - k 

d k=\ dtfc 


(15) 


— ^n— 1 [f] {.to, . •., t n , /A, z) 

where the conditions f n< q — 0 are implied. The solution of this equation should be sought in the 
following form 


f„=fn(t 2 ,t 3,...;Ai) +/n(fO,fl,---;A0 (16) 

where f n and f n are chosen such to satisfy, respectively, the following two equations: 

H - d-D (n) (1 - K 2 ) H [/) «o,. ■ •, z) (17) 

and 

- (jU) (1 - JU 2 ) = 2n- 1 [/] (f 2 , • • • T»; JU, z) (18) 

The solution for f n is as follows 


fn=Z c l n) (to)e-** to ilf k (n) (19) 

k=\ 

and it can be evaluated for arbitrary n by expanding both sides of eq.dTTI) in a series of 
eigenfunctions of the diffusion operator on its l.h.s.: 

For D — 1, for example, 1 / 4 . are the Legendre polynomials with — k(k + \ ), k — 0,1,_The 

in) 

time dependent coefficients C k are determined by the initial values of f„ (anisotropic part of the 
initial CR distribution) and the r.h.s. of eq.tfTTT). that depends on /„_ 1 , obtained at the preceding 
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step. It is seen, however, that f n exponentially decay in time for t > 1 and we may ignore therrQ 
as we are primarily interested in evolving the system over times t > e 2 1 and even longer. 

Starting from n — 0 and using eq. (fl5T) . for the slowly varying part of / we have 

§|=°- ( 20 ) 

The solubility condition for f\ (obtained by integrating both sides of eq. lfl5l in /i) also gives a 
trivial result 


dfo 

dt\ 


= 0 , 


( 21 ) 


so the last two conditions are consistent with the suggested decomposition in eq.(fl6l)- since from 
eq.([l8l) with n — 1 we have 


h = ~\ W ^ ( 22 ) 

and, thus both /o and f\ are, indeed, independent of to and t\ . We have introduced the function 
W (/i) here by the following two relations 


dW 

dfl 




(23) 


The solubility condition for /? yields the nontrivial and well-known (e.g., Jokipii 196 6 ) result, 
which is actually the leading term of the dfo/dt expansion in e <C 1 


where 


dfo 

dt2 


1 

4 




(24) 


K 



\ lie solubility conditions for fo, / 4 ,... will generate the higher order terms of our expansion 
which, after some algebra, can be manipulated into the following expressions for the third and 
fourth orders of approximation 


! In fact we must do so, as our asymptotic method has a power accuracy in e 1, but not the 
exponential accuracy. 
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dfo 

dt-i 

dfo 

dt 4 



1 

8 \dz 


+ 



x 




[k( 1-^) + U] 2 \\ df 0 
D( 1-M 2 ) // dz' 


(25) 


(26) 


We have denoted 


U = 


IzE 

D 


d\l, 


and U' — dU/dfi. The pitch-angle diffusion coefficient D(\i) and magnetic focusing a are 
considered z- independent for simplicity, a limitation that can be easily relaxed by re-arranging 
the operators containing d/dz in eq.(l26l ). We can proceed to higher orders of approximation ad 
infinitum since terms containing ((l — /i 2 ) dj n /dfl) can be expressed through /„_|. /„_ 2 ,.... 
According to eqs. (l20H2rT) . of interest is the evolution of fo on the time scales 1 2 > 1 or t > £ 2 
so, as we already mentioned, the contributions of /„ (/i) to all the solubility conditions, similar 
to those given by eqs. (!24H26l) . have to be dropped (as they become exponentially small) and 
only f n (/ 1 )- contributions should be retained. Using eqs. (l20H21l24ll26l) to form the combinations 
£ n d n fo/dt n and summing up both sides, on the l.h.s. of the resulting equation we simply obtain 
dfo/dt (see eq. lTHll !. Therefore, the evolution of /o up to the fourth order in £ takes the following 
form 


dfo £ 
dt 


2 

= T* 


■Z2 r 


K - (MVk 2 ) - — [Ki {d") -K 2 d' z d z 


\dh 

J dz 


where d' z = d z + o, d" — d z + o/2, and 


(27) 


K, = (W>(k-U')), = (28) 

The above algorithm allows one to obtain the master equation to arbitrary order in £. By 
construction, in no order of approximation will higher time derivatives emerge, as has been 
devised by Chapman and Enskog. We have truncated this process at the fourth order, £ 4 . As we 
show in the next section, this is the lowest order required to relate the above result to the telegraph 
equation. It also gives the first non-vanishing correction to the standard CR diffusion model in an 
important case (/iW 2 ) = 0, which is fulfilled, in particular, for = D (/ 1 ). Higher order 

terms can be calculated at the expense of a more involved algebra, but we argue below that such 
calculations would not change the results significantly. 















12 


4. Comparison with Earlier Results. Recovering Telegraph Term 


In contrast to the tele graph equation given by eqs.(|4H5]), that has been derived by 
Litvinenko & Schlickeiserl ( 20131) using direct iteration of eq.© with no explicit ordering of 
the emerging terms, eq. (l27l) is derived to the £ 4 - order of approximation with an e n factor 
labeling each term. Yet, it has no second order time derivative which is inconsistent with eq.©. 
Below we demonstrate that eq. (l27l) can still be converted to the telegraph form, however, with 
additional terms absent from eq.©©. Although eq.(l27l) is obtained by a broadly a applicable 
Chapman-Enskog method, its reduction to the telegraph form below is more restrictive and should 
be taken with a grain of salt for the reasons we discuss later. 

Several versions of telegraph equation have been obtained using different methods 
but, unfortunately, many of them do not offer clear orderin g, as e.g., eg.© derived by 


Litvinenko & Schlickeisen 120131 In an earlier treatment by (lEarll 1 1973h . an eigenfunction 


expansion was truncated with no transparent assessment of discarded terms. As we mentioned 
already, many treatments do not systematically eliminate short time scales which are irrelevant to 
the long-time evolution of the isotropic part of the CR distribution. In principle, this is acceptable 
if the reduction scheme is based on an exact solution of the original equation, to include all 
required orders of approximation into the master equation. Such approach, along with a nearly 
exhaustive analysis of the previous work has been presented in ( Schwadron & Gombosi 19941) . 
Their treatment, however, is by necessity limited to a relatively simple D (ju) (i.e., power-law in 
H). 


To clarify the role of the higher order terms in eq. (l27l) . we note that the r.h.s. of this 
equation represents just the first three non-vanishing contributions from an infinite asymptotic 
series (in £ 1) which we would obtain by continuing the reduction process described in the 

preceding section. This series may or may not converge to some linear (integral) operator in 
z. From a practical standpoint, the maximum order term that needs to be retained is either the 
first non-vanishing term, or else it introduces a new property to the solution, such as symmetry 
breaking. Precisely the last aspect has been highlighted by Litvinenko & Schlickeisei ( 2013 ) who 
used the telegraph term to calculate the skewness of a CR pulse. From this angle, we examine the 
third and the fourth order below separately. 


Third order equation. To this order eq.(l27l) rewrites 


dfo = ydfo 


F K\ 


d 2 fo 


llW 2 ) 


d 3 fo 


(29) 


dx dz ' 1 dz 2 ' / dz 3 

We have introduced the slow time x = e 2 t/ 4, as a natural time scale for the reduced system, and 
the following notation 


V = a 




(tiw 2 ). 



























13 


Note, that using x instead of 1 makes the terms looking lower by two orders in £, but describing 
them in the sequel we will use their original order, as it stands in eq. dTTT) rather than eqs. (l30l) or 

m. 

Eq.(l29l) can be solved using a Fourier transform and integral representations of Airy 
functions. We consider some basic properties of this solution using the moments of Jo. In 
particular, it is seen from this equation that the skewness of a CR pulse, propagating at the bulk 
speed V, arises in this order of approximation. Indeed, upon a Galilean transform to the reference 
frame moving with the speed —V, z —> z' — z + Vx, the above equation rewrites 


dfo 

dx 


Kl 


d 2 fo 

dz' 2 


£ {flW 2 ) 


d 3 fo 

dz' 3 


(30) 


so the last term generates an antisymmetric component of fo {z ), even if fo (z') is an even function 
of z' initially. By normalizing /o to unity 


/o = 



—oo 


1 , 


and assuming the coefficients in eq.(l30l) to be constant, for the moments of /o (z',l ) 


^ = / ^ )Z ' nd7 '"' 

— oo 

we obtain 

—z' = 0, —z' 2 — 2ki, — z' 3 = 6kiz' + 6e {fiW 2 ). 
dx dx dx w / 

With no loss of generality we may set z' = 0 and, in addition, z' 3 — 0 at T = 0, so that the skewness 
of the CR distribution changes in time as follows 


S = 



6 e(juW 2 ) 
(z'o +2 k\x^ 


X 

~J]2 


where ’0’ at z! 2 refers to its value at x — 0. The skewness remains small and its maximum 


2 £(/xW 2 ) 
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is achieved at x — z! 2 / Kj. Not surprisingly, the skewness increases with decreasing initial spatial 
dispersion of CRs. Indeed, according to eq.©), a narrow f(z) generates strong pitch-angle 
anisotropy which, in combination with asymmetric pitch-angle scattering ((/ 1 VF 2 ) f 0), generates 
the spatial skewness of the CR pulse. For not too small x, the explicit form of the solution of 
eq.(l30l) may be easily written down by using, e.g., a Fourier transform in z' and the steepest 
descent estimate of its inversion 


C 


f{z, t) ~ -^=exp 


,/2 


4k x x 


1 - 


£Z 

2k\x 




(31) 


It is quite possible, however, that even this small effect does not occur because of the pitch angle 
scattering symmetry, that is (/ilk 2 ) = 0. In this case the solution remains diffusive and, to obtain 
corrections to it and to see where the telegraph term might come from, the next approximation 
needs to be considered. 


Fourth order equation. The telegraph term. In the absence of e 3 terms, that is when 
(l~lW 2 ) = 0, eq.(l27l) takes the following form 


dfo 

dx 


+ £ 2 (K l -K 2 ) 




d 2 fo , w dfo 

i<2^rr + V2^r- + 


dz 2 


dz 


ff(£ 4 ) 


(32) 


o 2 £ 2 f 5 

*=*-—{ 4 *-* 

This equation, obtained within Chapman-Enskog method, and the telegraph equation ©, obtained 
by a direct iteration method, differ from each other in the second term on the l.h.s. The remaining 
terms of the two equations are equivalent even though not identical due to the insignificant ~ e 4 
corrections included in the coefficients k 2 and V on the r.h.s. of eq. (l32l) . 

To understand how the conflicting terms on the l.h.s. of both equations are related, we note 
that within the regular ordering scheme leading to eq. (l32l) the term in question must remain 
small, being nominally an £ 4 -term. However, as the conflicting terms in both equations are the 
higher-order derivatives, they may stick out from their order of approximation if the solution 
strongly varies in space and time. In order to preserve the overall solution integrity in such events 
a multi- (time)-scale or matched asymptotic expansion method is normally applied. We will argue 
that the telegraph equation approach to the CR transport does not handle this situation properly, 
as opposed to the Chapman-Enskog approach. But this is not to say that the two terms cannot be 
mapped to each other, when they are well in the validity range in the above sense. Note that in a 
number of other treatments the telegraph term was tacitly handled as one of the dominant terms. 
We pointed out that indeed, the term in question on the l.h.s of eq.(l32l) is small only insofar as 
the z- derivative does not change its order of approximation due to strong inhomogeneity, whose 


, V 2 = a 


K 


0 2 £ 2 


-K X 
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scale should not be less than the m.f.p., A. So, assuming strong inequality A- 1 (or simply 
d z ~ 1 in our dimensionless variables), we can express the high order spatial derivatives using a 
zero order (£ —* 0) version of this equation with sufficient accuracy. The result reads: 


dfo £ 2 
dz 2 k 2 


(Ki-K 2 ) 


d 2 fo 

dT 2 



G 2 £ 2 


8 



d 2 fo 
dz 2 


+v 2 


dfo 

dz 


+ G (e 4 ), 


(33) 


This equation is indeed equivalent to the telegraph equation by its form, but the equivalence 
requires not only £ <C 1 but also smooth variation in z and T, as not to raise the actual value 
of these terms significantly. Under these equivalence conditions, both the telegraph and the 
hyperdiffusive transport terms are just the corrections and may be safely ignored (especially if 
the £ 3 contribution is not empty). On the contrary, when the higher derivatives strongly enhance 
these terms, the equations cannot be mapped to each other and their solutions are disparate. One 
of them (or even both) may become less accurate than the underlying leading order (diffusive) 
approximation. This is quite common situation in asymptotic expansions when the form of the 
next order term should be selected on the ground of the least possible singularity it introduces 
into the expansion (cf. small denominator, secular growth etc. in mechanical problems, where the 
higher order approximations, if handled blithely, only aggravate disagreement with true solutions). 
The telegraph term correction to the diffusive approximation appears to come from that variety, 
as it generates singular components (5- and Heaviside functions) that are not only inconsistent 
with the strong pitch angle scattering and resulting spatial diffusion, but with the scatter-free limit 
of the parent differential equation itself. Therefore, the singular part of the telegraph solution is 
inherited from the derivation of telegraph equation. We compare the telegraph and hyperdiffusion 
type corrections to the basic diffusive propagation somewhat further in the next section. 


5. Relation between Telegraph and Hyperdiffusion Approximation 

We start with a relatively minor aspect of the differences between the two models. As 
we stated in SccjT] the telegraph coefficient in eq. (l33l) is inconsistent with some of the earlier 
derivations. In the simplest case D — 1, for example, after proper rescaling of T and z, it 



in the general case of D (—/i) f D (ji ), the appropriate equation for describing CR transport is 
that given by the lower, £ 3 - order, not included in eq.©. 


More importantly, the telegraph version of eq.(l32l) given by eq.(l33l) is valid only if the 
telegraph term (<* £ 2 , fourth order term) remains small compared to the other terms and the 
original ordering in eq.(l32l) is not violated by strong variations of the solution in space and time, 
as we pointed out earlier. We signify this by the “slow” time z ~ e 2 t. In most other treatments 
t is used instead, which formally makes the telegraph term in eq.(l33l) appearing as a zero order 
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term. It is not important, of course, whether the term is labeled by e 2 or not; important is that it is 
treated as a subordinate term. Attempts to make it dominant a posteriori violates assumptions that 
are essential for its derivation. This point is demonstrated below by repeating a simple calculation 
of the CR pulse skewness that we already made earlier working to e 3 order. 


Litvinenko & Schlickeiserl ( 2013b suggested to study an asymmetry (skewness) of a CR pulse 


propagating along the field under the action of magnetic focusing using the telegraph equation. 
This and some other characteristics of the CR pulse, such as the kurtosis, can be easily analyzed 
using the primary equation (l32l) . The calculation of the pulse skewness essentially repeats the one 
already done at the e 3 level, where it is generated by asymmetric scattering, D ( —ji) / D (p). So, 
transforming eq.(l32l) to the reference frame moving with the speed —V 2 , that is z —* z! — z + V 2 T, 
we obtain 


^-z' = 0, ^- z ' 2 = 2k 2 , ^-z' 3 = 6£ 2 o(K 1 -K 2 ), (34) 

a T a T a % 

where we have, again, assumed z! — 0. The skewness thus evolves in time as follows 

c 6e 2 (j -K 2 ) t 

b \ 3/2 

\Zq + 2K 2 ZJ 

Unless S(z) reaches its maximum very early it is fairly small. Because T max = Zq/k 2 , the 
maximum value 


3'max — S (T m ax) — 


2e 1 o( K K\ —K 2 ) 


k 2 


3 Zn 


(35) 


so that an initially symmetric CR pulse develops significant asymmetry only if z!q < £ 4 . This, 


however, would require T max ~ £ 4 (or, equivalently, / n 


£ 2 ), in strong violation of the 


requirement t > 1, established in Secj3] Note that significant pulse asymmetry obtained by 
Litvinenko & Schlickeiserl ( 201 3| ) using the telegraph equation was based on the fundamental 
solution to this equation, that is z' {) —> 0. We argued in Secfflthat the early propagation phase is 
not adequately described by the telegraph equation so that the pulse asymmetry might have been 
overestimated in the above paper. We specify the validity range of the telegraph equation below. 

Starting from a long time regime z > £z r , similarly to the £ 3 result given in eq. (l3TT) . from 
eq.(l32l) we find 


C 


f{z,z) ~ y^e x P 


_/2 


4k 2 z 


, „ 2 (Ki-K 2 )z' 

2k 2 z 



(36) 
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This propagation regime is not much different from the regular diffusion (first term in the square 
bracket), so both the hyperdiffusion and telegraph models produce similar results, as they are 
largely equivalent in this regime. It is worthwhile to write the requirement for the agreement 
between the two models in physical units which is simply 


vt > z 

The point z — vt is close to the cut-off in the telegraph solution, z/t — \fkjf (for a = 0), eq.©. 
It follows then that unless z 2 3> kt (t T), the cut-off strongly changes the overall solution. 


The opposite case T < £z', which corresponds to the initial phase of pulse relaxation, is the 
key to understanding the difference between the Chapman-Enskog and the telegraph methods. 

In this regime they are not equivalent , as the e 4 term in eq.(l32l) cannot be neglected to make the 
transition to the telegraph equation (l33l) . Indeed, when propagation starts with an infinitely narrow 
pulse, in the early phase of its relaxation the higher z- derivatives are still too large for such 
transformation. A spatially narrow pulse automatically generates strong pitch angle anisotropy 
which, in turn, results in rapid time variation, making the telegraph term also large. It is this 
regime where both methods become questionable and, in addition, their predictions deviate from 
each other both quantitatively and qualitatively. We need to check first whether they are relevant 
to this regime. 


The phase t < £z! is well described by the Chapman-Enskog approach down to z ~ 1/e 2 
(to ~ 1, Sec©. The situation with the telegraph equation is more complex, as the conversion from 
the Chapman-Enskog expansion is invalid, while independent derivations rarely provide clear 
ordering. A rigorous derivation in ( Schwadron & Gombosi 19941) requires the same assumptions 
that we made when transforming the hyperdiffusive equation into the telegraph equation, that is 
c) q- rvj d 2 (e t rs./ £ 2 under their nomenclature). So, the telegraph equation appears to be a subset 
of the hyperdiffusion equation valid only under the above ordering. It is likely to break down in 
the r < £z! regime, in other words, near its cut-off. We support this premise by the following 
considerations. 


Recently, Effenberger & Litvinenko ( 20141) and Litvinenko et al. ( 2015 ) have carried out 
simulations of the full scattering problem, corresponding to eq.®. The results deviate from the 
telegraph solution precisely at the early phase of the pulse propagation, when the hyperdiffusion 
and telegraph models disagree. The two 8- function pulses with sharp fronts in its solution given 
by eq.© are not seen in the simulations. This is understandable, as such features are inconsistent 
with the underlying scattering problem. They should have been smeared out by scattering earlier, 
since the spatial profile is shown at five collision times (Fig.2 in Effenberger & Litvinenko! 
2014). Moreover, the 8— function pulses 8 ( z ± \/k/T t j that are an integral part of the telegraph 


solution, as they maintain its normalization, are irrelevant to the primary equation ® , even without 
collisions. Indeed, if S> = 0, and the initial condition is 8 (z) being constant in /i for — 1 < /i < 1, 
the scatter-free solution is 8(z — fxvt). Hence, fo(z,t) = / fdfl/2 — (2vt)~ l H (vt — \z\). 
Therefore, a certain property of the telegraph equation allows the 8 (z ± \JkjT t ) and sharp 
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front components (although with modified speed, eq.[[6]]) to survive multiple collisions. As these 
components are inconsistent with the underlying scattering problem, this property of the equation 
must have been acquired during its derivation. Obviously, it is rooted in the hyperbolic (telegraph) 
operator df - (k/T) dr that allows the singular profiles to propagate without spreading. 

By contrast, the hyperdiffusion equation does not require singular components but, on 
the contrary, smears them out. We present an approximate solution of eq. (l32l) after neglecting 
magnetic focusing in e 4 - order terms in eq.(l32l). We also neglect the regular diffusion compared to 
the hyperdiffusion, which is acceptable during an early phase of pulse relaxation, z < £z. Finally, 
we assume z! > 0, as the solution is an even function of z'■ The asymptotic result is as follows: 




■1/6.' !/ 3 exp 


3 1/3 A 3 \ /3 3 / 2 

-g 4 ww ,cos 


4 - 1/3 


J4/3 


h l/3 T l/3 


(37) 


We have denoted the hyperdiffusion constant h = e 2 {K\ — Ki) /2. More about this result and 
further discussion of the two conflicting approaches can be found in Appendix. We see that 
there is a considerable slow down of the CR spreading compared to the conventional diffusion, 
z' 2 r 1 / 2 , that embodies a sub-diffusive propagation, z' 2 f l,/4 . This ameliorates the problem 

of acasual propagation in diffusion regime, yet no sharp fronts or spikes develop. By contrast, 
the telegraph solution to the causality problem is to cut off the solution beyond certain distance 
(|z| > y/k/Tt), thus introducing an unphysical singularity. The immediately arising normalization 
problem is then “solved” by adding an even stronger singularity in form of two 5 functions at the 
cut-off points. 


6. Summary and Conclusions 

Using the Chapman-Enskog method, we have extended the CR transport equation with 
magnetic focusing to the fourth order in a small parameter e = A // (CR mean free path to the 
characteristic scale of the problem). This analysis clarifies the nature of the telegraph transport 
equation, widely publicized in the literature as a promising alternative to diffusive propagation 
models. We have shown that the telegraph extension (°c d 2 fo/dt 2 ) of the diffusion equation can 
be mapped from the (small) hyper-diffusive term (°c <9 4 / 0 / dz 4 ) of the regular Chapman-Enskog 
expansion, but the telegraph term, originating from an e 4 term of the expansion, must remain 
subordinate to the main, diffusive transport and magnetic focusing (if present) contributions. This 
condition is met after vt > z/A collision times for an initially narrow (Az < A) CR distribution. 
Another important limitation of the telegraph equation is that, by contrast to the Chapman-Enskog 
equation and the original pitch-angle scattering equation, it is not self-contained requiring an 
initial condition for also dfo/dt, which needs information about the anisotropic part of the initial 
CR distribution. Furthermore, an attempt to proceed to higher orders in £ introduces progressively 
shorter time scales associated with “ghost” terms reflecting quick relaxation of initial anisotropy 
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or strong spatial inhomogeneity. By contrast, the classic Chapman-Enskog method is devised to 
eliminate short time scales, irrelevant to the evolution of the isotropic part of CR distribution /o, 
which accurately describes this evolution after a few collision times, vt > 1. 

We have derived the CR transport equation for an arbitrary pitch-angle scattering coefficient 
D (ju). This form of transport equation, (1271) . is suitable for describing CR acceleration and escape 
problems wher e the phenomenon of self-confinement (D is a functional of f, D = D [/;/I,/]) is 
critical, e.g. ( Ptuskin et al.1 2008 : Maikov et al. 201 Obi 20131: Fujita et al. 2011 ). Accounting for 
magnetic focusing effects is required, e.g., for describing particle acceleration in CR-modified 
shocks with an oblique magnetic field. In this case the field increases towards the shocks due 
to the pressure exerted by the accelerated CRs on the flow, thus producing a mirror effect. The 
particle drift velocity along the field associated with the mirror effect is (e.g., cq. If29l ) V r k/I b 
which, for the magnetic field variation scale being of the order of the shock precursor scale K/U s h 
and strong shock modification, almost automatically becomes comparable with the shock velocity 
U s h- This additional bulk motion of the accelerated CR (directed towards the shock) will affect 
their spectrum and acceleration time. 

Furthermore, as the CR scattering in such environments (i.e., supernova remnant [SNR] 
shocks) must be self-sustained by virtue of instabilities of the CR distribution (see, e.g., 

Bykov et al. ( 2013b : Belli ( 2014 ) for the recent reviews), the above magnetic drift needs to be 
included in the CR stability analysis. It should be noted, however, that the results obtained in 
the present paper formally require a magnetic field Bo that does not strongly change over the 
gyro-radius of energetic particles. This is not to be expected in SNR shocks, especially if strong, 
CR current- and pressure-driven instabilities generate fields with 8B > Bq. However, one may 
use the shock normal direction as the polar axis to calculate the pitch angle diffusion coefficient 
D(ju), needed for the description of the CR spatial transport. Also, such treatment will require a 
description of the gyro-motion and averaging by computing particle orbits beyond the standard 
quasi-linear description (Maikov & Diamondll2006h . implied throughout this paper. 


In conclusion, by comparison with the telegraph equation, the classic Chapman-Enskog 
hyper-diffusion equation consistently describes the long-term CR propagation in a self-contained, 
order controlled fashion. Further improvement of the CR diffusion models should probably 
address the anisotropic component of the CR distribution. There are situations, such as ultra-high 
energy CR propagation, where the mean free path grows too long with th e en ergy as to make 
the diffusive approach irrelevant and a rectilinear transport to dominate (Aloisio et al. 2009b (cf. 
Fevi flight regime described in aforementioned study by Maikov & Diamond 20061). Another 
interesting example is a sharp angula r anisotropy ~ 10° in CR arrival directions discovered by 


MIFAGRO observatory (Abdo et al. 2008) and a number of othe r instruments, e.g. (lAbbasi et al. 


201 ll: iBartoli et al.l 120131: iDesiatil 120 1 4l: lAbev sekara et al.ll2014l). Unless this anisotropy is 


of a very local origin (such as heliosphere, dFazarian & Desiatil 120 1 Ol: lO’C. Druryll2013b . it 
poses a real challenge to CR propagation mo dels and clearly cannot be addressed within the 
diffusive approaches discussed in t his paper dDrurv & Aharonianll2008l: iMalkov et al.ll2010al: 
Giacinti & Siglll2012l: Ahlersll20 14 : Malkovll20 15 ). On the other hand, when the diffusive transport 
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model is well within its validity range (weakly anisotropic spatially smooth CR distributions) 
neither the telegraph nor the hyper-diffusive term (both~ £ 4 ) is essential to the CR transport and 
can be neglected. 


Useful and stimulating comments of the Referee are gratefully acknowledged. The work 
of MM was supported by the NASA ATP-program under the Grant NNX14AH36G. He is also 
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A. Appendix 

To derive the result given by eq.(l37l). we rewrite eq.(l32l) for a simple case of weak focusing 
(of < d z f) and using the Galilean transform to the frame moving at the magnetic drift velocity 
V 2 ,z' = z + V 2 z: 


df = K d2f l dAf 

dz 2 dz' 2 dz' 4 ' 


(Al) 


where h = £ 2 (K\ — K 2 ) /2 and k 2 is the same as in eq. (l32l) . The fundamental solution of cq. dA 1 1 ) 
can be written as an inversion of the Fourier image /*. (t), assuming that (0) = \/2n. As we are 
going to find the Fourier inversion using asymptotic methods, we write an arbitrary constant C 
instead of this value: 


f(z,z')=C J e ikzl ~ k2 ( K2+hk2 ) T dk. (A2) 

—oo 

We will determine the normalization constant C shortly from the requirement / —>■ 8 (z!) , T —> 0. 
It is convenient to introduce the following notation: 

£ 4 = 4hzk 4 , ^ = z! (4hz)^ 4 > 0. 

We may limit our consideration to the ^ > 0 half-space, because the solution is an even function 
of z!■ Focusing on a short time asymptotic regime z < £z!, which is opposite to the case considered 
earlier (eq. l[36l f we neglect the diffusive term ~ k 2 in the exponent (as hyperdiffusion dominates) 
and rewrite eq. (IA2l) as follows 


oo 
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Note that the general evaluation of the integral in eq. (IA2l) . with the diffusive term included, is not 
difficult but more cumbersome. The phase of the integral here has three saddle points, and the 
following two should be on the integration path 


f ± =i^/3 e ±«/3, 

since the integrand reaches its maxima at these points. So, the integration runs from — °° through 
£ + to +/°° then to £_ and, finally to +°°. The contributions from the two saddle points then yield 

/= i^ exp (-^ 4/3 ) cos (^ 4/3 ) <A4) 

The normalization f fdz! — 1 requires the constant C — 1 / V3k, which only insignificantly 
deviates from the exact value C — \/2n that, in turn, follows from the integral representation of 
/ in eq. (IA2l) for T —> 0. Note that we could have replaced an oscillating exponential tail of this 
solution by zero beyond the point E, > (47r) 3//4 3 -9 / 8 . Such modification of the hyperdiffusive 
solution would be in the spirit of the telegraph cut-off at z/t — y/k/T, however, with an essential 
difference of being only a discontinuity in the solution derivative. The unphysical (oscillatory) 
behavior beyond the first zero point of the solution in eq. (lA4l) results from neglecting the diffusion 
term in eq. dAll) . asymptotic methods used to calculate the integral in eq. (lA3l) . and from lacking 
higher order terms, ~ e'\ n > 4. Therefore, the solution can be improved systematically. In 
addition, it starts from a point source which is clearly inconsistent with the main approximation 
£C 1. A somewhat broader initial profile will not develop an oscillatory tail, if convolved with 
the Green’s function in eq. (lA2l) . We will not attempt to improve on this minor aspect of the 
solution here, as it becomes only weakly irregular if cut off at its first zero. 

From the perspective of a general improvement of the asymptotic expansion considered in 
this paper, the derivation of eqs.(l34l). for example, is robust in the following sense. The residual 
higher order terms in e -C 1, if included in eq.(l32l). will not change eqs.(l34l) in any other way 
than small corrections to the coefficients K2 and Vi- Indeed, the higher z- derivatives coming 
from higher order terms, will vanish from the (first four) moment equations after integrating by 
parts. By contrast, continuing the telegraph approach to higher orders will generate terms with 
small parameters at higher time derivatives in all moment equations. These terms will become 
crucial during the initial relaxation of the CR distribution. The relaxation is associated with the 
CR anisotropy or strong initial inhomogeneity, that is with large /„, Secj3] However, these decay 
over a short time t < 1 . This is the time period when the telegraph or hyperdiffusive correction 
is large but its effect on the subsequent evolution ought to be limited as this time is short. The 
hyperdiffusive correction meets this requirement, as we argued using moment equations. To see 
whether the same is true for the telegraph correction, let us rewrite eq. (l33l) using the “fast” time, 
t = t / c 2 : 
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+ T dt) dt 4 dz 2 '' 


(A5) 


where we denoted Tt = 2 (K\ — Kf) /k ~ 1 and assumed <7 = 0, to make the following simple 
argument. Namely, in the limit £ —y 0 there are two modes, of which the first is fo = fo (z). This is 
the main diffusion mode that slowly evolves in time when 0 < £ -C 1, and, as we are interested in 
the evolution over the time scales t > £ 2 . the telegraph term becomes ~ £ 4 and can be discarded. 
The second mode corresponds to a rapid decay of the initial distribution ~ exp (— t/xj) which 
is associated with the decay of initial anisotropy or strong inhomogeneity. If this mode is active 
(dfo/dt^ 0 in eq. llA51 1. then even the total number of particles N is not conserved automatically. 
So, turning to the moments of eq. (lA5l) we need to impose the initial condition, dN/dt — 0, 
to ensure the particle conservation. This probably means that £ —» 0 is a difficult limit for the 
telegraph reduction scheme. The singular components in the telegraph solution © appear to be 
primarily associated with the particle conservation problem. The initial relaxation phase {t < 1) 
perhaps, cannot be adequately described by the telegraph reduction scheme using an equation 
for fo alone, as it does not properly “average out” an anisotropic component /, which is large 
during this period of time. The telegraph term is therefore to be understood as a “ghost” term 
reflecting rapid decay of such components. It follows that the rapidly changing part / in the 
decomposition in cq.dTbl) needs to be retained in the short-time analysis along with fo. Otherwise, 
the telegraph operator generates unphysical 8- pulses and sharp fronts, just to conserve the 
number of particles, as discussed earlier. These considerations are, however, not nearly complete. 
Further useful analysis of propagation modes in the context of telegraph equation can be found in 
( Schwadron & Gombosil 1994 1. 


To conclude this Appendix we make yet another argument in disfavor of the telegraph 
equation that is partially related to the above considerations. A consistent asymptotic reduction 
method must be continuable to infinity in powers of small £. The Chapman-Enskog scheme 
clearly is. The outcome will be a series of terms ~ d"/o on the r.h.s. of eq.(f27l) with just d t fo 
on its l.h.s. To solve the resulting equation, only the initial distribution fo (0,z) is needed, as the 
equation remains evolutionary and (generalized) parabolic, as its pitch-angle diffusion superset 
is. The telegraph equation, on the contrary, turns hyperbolic and non-evolutionary after the 
reduction from the superset equation. By continuing the telegraph reduction scheme to higher 
orders of approximation, progressively higher time derivatives will emerge (along with higher 
space derivatives). The resulting equations will thus be non-evolutionary and a growing set of 
initial time derivatives d n fo will then be needed to solve the initial value problem. These data 
can be extracted only from the full anisotropic distribution with recourse to the full (anisotropic) 
equation. Therefore, the telegraph equation is not self-contained and cannot be improved 
systematically. Any attempts to improve it will introduce shorter and shorter time scales that 
would require a return to the full anisotropic description. 
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